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We review the phenomenology of vr" production at large transverse momentum 
0^ ' in proton- induced collisions. Uncertainties in the next-to-leading-order predictions of 

^^ I Quantum Chromodynamics are discussed. The comparison with data reveals that the 

disagreement between theory and experiment lies essentially in an overall normaliza- 
C^h! tion factor. The situation for tt'^ production is contrasted with that of prompt-photon 

(->! • production in hadronic collisions. 
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1 Introduction 

The production of hadrons at large transverse momentum pT (in hadron-hadron or hadron- 
photon collisions) offers a classical test of perturbative QCD. The cross sections of the 
hard subprocesses are calculated with next-to-leading-order (NLO) accuracy |||], the quark 
and gluon distributions in the initial hadrons are measured in deep-inelastic-scattering 
experiments |^,y,^j5| and the fragmentation functions describing the transitions of the 
partons into the final-state hadrons are extracted from e^e~ -annihilation data ||6|,0Q. 
Therefore, all the building blocks of the large-pT cross sections are in principle known, 
and the comparison with data provides interesting tests of the theory. 

The hadroproduction of large-pr vr'' mesons has been studied by the CGGRW Collab- 
oration [^, who obtained good agreement between the QCD predictions and collider vr*^ 
cross sections. On the other hand, fixed-target data with a lower center-of-mass energy y^ 
(23 GeV^y^^SO GeV) overshoot the theoretical predictions, and it was impossible to 
obtain an overall agreement with all the large-px vr*^ data. Lavge-pT charged-hadron cross 
sections were studied in hadroproduction at collider energies Q and in photoproduction 
at HERA energies ||^]. In that case, reasonable agreement between theory and data was 
reached, but one must remember that these predictions are sensitive to the choice of the 
factorization and renormalization scales in the pr-i'ange which was studied {px ^ 10 GeV). 

In this paper, we come back to the study of the laige-pT vr" cross sections. New fixed- 
target data appeared recently |10,11[, completing those already available in the range 



23^ V^^30 GeV |T^. On the theoretical side, new sets of fragmentation function ex- 
tracted from LEP data are now available [Ml . One must also notice that A^-J- has increased 

LLH MS 

from ~ 200 MeV to ~ 300 MeV since the time of the first study ^] . This new value of Aj^ 
gives rise to a non-negligible increase of the QCD cross section in the pT range studied 
here {pr^W GeV). 

Another important reason for this study is the publication of new data on prompt- 
photon cross sections [^,^. A recent QCD analysis of all fixed-target and ISR data [13| 



leads to the conclusion that theory is in reasonable agreement with experiment, with the 
exception of two data sets, as shown in Fig. 0. However, it should be remembered that 
the comparison with QCD predictions can only be done for high-enough values of px (s-g- 
Pt^5 GeV/c for the E706 and ISR energy ranges), where the scale dependence of the 
theoretical cross sections is not too large. Therefore, in Fig. |l|, one should not consider the 



first points of the E706 and ISR experiments. Tlie introduction of an additional intrinsic 
transverse momentum kx of the incoming partons, strongly advocated in Refs. pO|JI^ to 
enforce agreement between QCD predictions and experiment, is not a very satisfactory 
solution; the agreement with the E706 data sets is improved, but the agreement with other 
data sets is reduced. 

This situation is challenging and leads us to look at large-pT '''"'' production. The 
large-pT t^^ mesons form a significant background for the prompt photons, and their 
cross section must be carefully measured in order to allow for a reliable estimate of the 
"fake prompt photons," in particular due to configurations in which one photon from 
the decaying vr'^ escapes detection. Therefore, there is a strong experimental correlation 
between the prompt-photon and the t:^ cross sections, especially at low px values, where 
the vr'^ background is largest. A study of the latter could bring some clarification on the 
prompt-photon puzzle. 

Section 2 of this paper is devoted to a discussion of the theoretical calculation of the 
vr'^ cross sections. We emphasize the uncertainties associated with the QCD scales, the 
determination of the fragmentation functions and the importance of the higher-order (HO) 
corrections. We discuss the comparison of theoretical predictions with data at fixed-target 
energies in section 3 and with data at higher energies in section 4. Section 5 contains the 
conclusions. 

2 Theoretical Framework 

At NLO, the inclusive cross section for the hadroproduction of a single hadron /i, differ- 
ential in the transverse momentum px and the pseudorapidity r] of h, takes the following 
form: 



i,j,k=q,g 
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The parton densities of the incoming hadrons hi and /12 are given by -Fj//ji and Fj/f^^ ; 
the fragmentation of a parton k into the hadron h is described by the fragmentation 
function D^{z, Mp); daij^k/dpTdrj is the Born cross section of the subprocess i+j — > k+X; 
and Kij^k is the corresponding HO correction term. In this paper, we shall use the ABFOW 



[^, CTEQ |3|,^] and MRST [^ parton densities. For the fragmentation functions, we shall 
use the paranietrization of Ref. Q. These fragmentation functions were derived from 
fits to charged-pion spectra, but we make the usual assumption, supported by data, that 
the rate of vr'^ production is half of that for charged pions. They are preferable to the 
parametrizations of Ref. [^, which predate the publication of the LEP hadronic spectra. 
Expression (ffl) depends on the initial- and final-state factorization scales, M and Mp, 
and on the renormalization scale /x. The cross section da^/dpTdrj, calculated to all orders 
in as, is independent of M, Mp and fx. However, the perturbative series calculated at 
fixed order in Os does depend on these scales, the compensation between the variations of 
Fijh^{M), Fj/fi,^{M), D^{Mp) and as(/i) and those of Kij^hil^^M^Mp) being incomplete. 
Therefore, we shall carefully explore the scale dependence of the cross section (|l]). On the 
one hand, we shall use "standard scales" M = Mp = /i = Pt/ K', with k, varying around 
K = 2, and, on the other hand, we shall fix the scales by using the Principle of Minimum 
Sensitivity (PMS) H, 



d_ da^ _ _d_ da'' _ d da'' 

djjL dprdr] dM dpTdt] dMp d-ppdr] 

which determines the optimal scales fj-opt, Mapt and Mp^opt- An illustration of the scale 
sensitivity is presented in Fig. ^, where the differential cross section for the inclusive pro- 
duction of single ir^ mesons, with px = 7 GeV/c and |r/| < 0.75, through the scattering of 
530 GeV/c protons on a fixed Beryllium target is considered. Specifically, the contours of 
constant cross section are shown in the {M,Mp) plane. For each set (M,Mp), /i is deter- 
mined by the first equality in Eq. (|2|). As expected, one finds a minimum, corresponding to 
the optimal choice defined by Eq. (|2|) , but one also notices that the area of stability under 
changes of scales around this optimal point is not very large. It is smaller than in the case 
of prompt-photon production. We noticed in several numerical studies that it is not pos- 
sible to verify the PMS criterion (^) for too small values oipp, typically for pr ^ 5 GeV/c. 
This indicates that the HO corrections are large and that the QCD predictions are quite 
unstable under changes of scales. Therefore, we shall emphasize the comparison between 
theory and experiment at pT^5 GeV/c, although we shall also explore smaller pp values. 
Since working in the optimal scheme is rather cumbersome and requires a lot of numer- 
ical work, we compare, in general, theory predictions for fixed scales with experiment. We 
use the standard choices M = Mp = fJ- = prf^ and pr/^- As will be discussed later, the 
optimal scales, derived from Eq. (pi), turn out to interpolate between these two choices. 



depending on the kinematical range considered. 

Besides the scale dependence, there are other sources of uncertainties in Eq. (|I|), which 
come from the parton densities and the fragmentation functions. The parton densities 
are generally determined with good accuracy. However, at small pT and for scale pt/'^, 
the factorization scale may approach the starting scale Qq of the QCD evolution. In this 
region, the parton densities are not constrained by data. To avoid this problem when 
using the scale pt/3, we thus require pT to be larger than 4.5 GeV/c. 

Another source of uncertainty concerns the fragmentation functions. The quark frag- 
mentation is now well constrained by PETRA, PEP and LEP data |0,3- However, the 
dominant support to Eq. (|^) comes from the large-z domain, where the e^e~ data are 
scarce. For instance, numerical studies of Eq. (||) indicate that the mean z value is about 
0.85 for 4 < pt- < 8 GeV/c and ^/s = 31 GeV. There are very few data points in the range 
0.8 < z < 1 (zero or one point in most data sets), and they have large error bars. Therefore, 
we expect the uncertainty coming from the fragmentation functions to be of the order of a 
few tens of percents. It must also be noticed that the gluon fragmentation function is not 
well determined by e^e~ data, since it appears there only at NLO. More constraints could, 
in principle, be obtained from inclusive pion production at large transverse momentum 
in hadronic collisions. For example, the inclusive large-p^ charged-hadron cross section 
has been measured by the UAl Collaboration [16| in the range 5^p7^^20 GeV/c. These 



data are compatible with the parametrization of Ref. H|, but they are also compatible with 
the same parametrization after the normalization of the gluon fragmentation function has 
been increased by 30% to 40%. (Such a large change in the gluon fragmentation function 
is still compatible with the e'^e~ data.) This gives us an estimate of the flexibility in the 
gluon parametrization. We remark that changing the gluon fragmentation function may 
affect the slope of the pT distributions, since the low xt regions will be more strongly af- 
fected. Besides these problems related to the z behavior of the fragmentation functions, we 
should remark that we are led to explore small fragmentation scales (in particular, when 
using Mp = pt/'^) fai" away from the kinematical regions where the fits are performed. In 
summary, hadronic data require large z values and small scales, while the fragmentation 
functions are (mostly) extracted from e~^e~ data at medium z values and large scales. 

Finally, let us mention that we do not use theoretical expression in which the large 
ln(l — z) terms present in the HO corrections are resummed ll^JT^]. Indeed, in order to 
be coherent, one should also use resummed expressions when extracting fragmentation 



functions from e'^e data. Such an analysis has not been done so far. It would certainly 
be very interesting to pursue phenonienological studies in this direction [19|. One must, 



however, keep in mind that the optimization procedure of Stevenson and Politzer [15| 
amounts to a partial resummation of the ln(l — z) terms [pO| . 

For all these reasons, we do not expect a very good agreement between fixed-target 
data and QCD predictions. The expected agreement should only be within a few tens of 
percent. At larger energies, the mean value of z is smaller and the sensitivity to the scale 
variations is reduced. A better agreement is therefore expected and has been verified in 
previous studies [y,y]. This point will be discussed later. 

In this paper, all the calculations are performed in the MS scheme. The value of Aj^ 
in Os, calculated at two loops, is taken to be equal to the one used in the parton densities. 

3 Fixed- Target Data and Theory 
3.1 Inclusive Distributions in Pt 

The data that we are going to discuss are displayed in Fig. y. Unlike what was observed 
in the case of prompt-photon production, the vr^ data taken at the same energy are quite 
compatible with each other, with the exception of ISR data at ^/s = 63 GeV, where the 
R806 [^] and AFS |22] data disagree for pt^6 GeV/c. The increase in cross sections as 



-^s increases at fixed pT is also visible. Let us emphasize that we only consider pp, pp and 
pBe data in order to reduce the uncertainties due to the parton densities of the incoming 
hadrons. We do not study irp data, but concentrate on proton-induced fixed-target data 
and ISR data in the energy range below ^/s = 63 GeV. 

In what follows, we consider only the range pr > 4 GeV/c, which is covered by 
almost all data sets. With this cut, we also efficiently suppress possible non-perturbative 
contributions to the cross sections, such as power corrections, intrinsic /ct effects, etc. 

In Fig. H, we show the effects of the variation in the parton densities and scales on 



the predictions for the WA70 experiment |12]. For both scale choices, pt/2 and pr/3. 



we observe the following hierarchy: ABFOW and CTEQ5M lead to the smallest and 
largest predictions, respectively, while MRST2 leads to an intermediate result. The result 
for CTEQ4M almost coincides with that for MRST2 and is not shown in Fig. |^. This 
hierarchy may be explained by the different Aj^ values used in the evaluations of Og. For 
four quark flavors, these are 230 MeV for ABFOW, 296 MeV for CTEQ4M, 300 MeV 



for MRST2 and 326 MeV for CTEQ5M. When Og is evaluated with a common value of 
Aj^, the four results become very similar. The effect of the scale choice is also clearly 
displayed. One notices a pT-independent increase by roughly a factor of two when the 
scale is reduced from pt/2 to pt/3. This is a typical example for the sensitivity of the 
theory to the scale choice. As we shall see later, the results obtained using optimal scales 
are close to those obtained with scale pt/3. As expected, the agreement between data and 
theory is not very good. With small scales (~ Pt/3), the theory underestimates the data 
by some 40% to 50%. In Fig. @, UA6 pp data lO] are compared with the corresponding 
QCD predictions. For this experiment, the QCD predictions for scale pr/3 undershoot 
the data by some 30%. 



Similar remarks hold for the E706 data [10|. In Fig. g, we display data and predic- 



tions for y/s = 31.6 GeV. One observes that the curve corresponding to the optimal scales 
(labelled PMS) is very close to the one obtained with scale pt/3. Here again, the dis- 
agreement between theory and data, for scale pt/3 or optimal scales, is of the order of a 
few tens of percents. It is interesting to notice that, with scale pt/2, which allows one to 
extend the predictions down to smaller values of pT, the theoretical cross section is almost 
parallel to the data. Therefore, there is no such increase of the ratio of data over theory 
as was observed in the prompt-photon case (see Fig. |l|). Very similar conclusions can be 
drawn from the comparison of the E706 data at ^/s = 38.8 GeV with theory. 

Finally, in Fig. U, we show the ISR R806 data at ^/s = 30.6 GeV [21| and the theoretical 



predictions obtained with scales pt/2 and pt/'^- Here again, small scales lead to the best 
agreement with the data. 

From this short survey of data and QCD predictions we can draw some preliminary 
conclusions. As expected, we do not obtain a very good agreement with respect to nor- 
malization. The predictions for scale pt/3 systematically underestimate the data by some 
30% to 50% for all considered experiments. On the other hand, the pT behavior is well 
reproduced. One also notices that all data appear to be consistent with each other. 

3.2 Inclusive Distributions in xt 

The above discussion is best summarized by displaying on a linear scale the data as differ- 
ential cross sections in xt normalized to the theoretical predictions. The data, normalized 
to the NLO predictions evaluated using set MRST2 and a common scale set equal to 
Pt/2, are displayed in Figs. ^ and ^. We notice that the normalized data are mutually 



compatible within ±20%, and that multiplying the theoretical predictions by a common 
normalization factor of 2.5 will bring data and theory in reasonable agreement. It is also 
interesting to remark that the ratio of data over theory is rather flat and shows no such 
sharp rise at low xt as was observed in the prompt-photon production experiments (see 
Fig. I). 

These general features are also found when comparing data and theory with the com- 
mon scale put equal to pt/3. The results are shown in Figs. 10 and O for the MRST2 



parton distributions. The major difference between theory and experiment now resides in 
a normalization factor of about 1.45. Comparing Fig. 11 with Fig. ^ one notices the differ- 



ent slope of the data-over-theory ratio at small xt values, corresponding to p-r ^ 5 GeV/c, 
indicating that the theoretical predictions are rather unstable in this region and cannot 
be trusted. 

3.3 Discussion 

Based on the above analysis of fixed-target data (WA70, UA6, E706, R806 at ^/s = 
30.6 GeV), we reach the following conclusions concerning inclusive single- vr*^ production. 
All data appear to be consistent with each other to within ±20%. However, they differ 
from the fixed-scale theoretical predictions by a rather large normalization factor, namely 
K = 2.5 ± 0.5 for scale pt/2 and K = 1.45 ± 0.25 for scale pt/'^- In the latter case, it is 
worth investigating if this normalization factor could be accounted for, in the theoretical 
calculations, by a different choice of fragmentation functions, which are still rather flexible 
in the kinematical range of interest, as explained in section |2|. 

This situation is to be compared to that of prompt-photon production, where the 
scattering of data sets is much larger. In fact, normalized to the theoretical predictions 
for fixed scales, the measured E706 rate at ^/s = 31.6 GeV appears two to four times 
larger than the corresponding WA70 rate. 

Another feature distinguishes the vr" spectra from the prompt-photon spectra when 
they are compared to theory. For scale pt/2, the data-over-theory ratio for the vr'' data is 
rather flat down to rather low xt values (see Fig. ^, while the corresponding ratio for the 
prompt-photon data is found to exhibit a sharp rise for decreasing xt- For scale pt/3, a rise 
is observed for both vr'^ and prompt-photon data, the rise being sharper in the latter case. 
This change of behavior signals an instability of the theory: for vr'^ production at E706 
energies, a change of scales introduces a change of slope in the pT spectrum. For WA70 



or UA6 energies, the change of scales reduces to a change in the overah normaHzation. 

A final remark concerns the energy variation of the n^ spectrum. By comparing data 
at different energies in the very same experiment, we are able to better examine the 
variation with ^/s of the theoretical predictions. This can be done by looking at Figs. |^ 



and 11, where E706 data at ^/s = 31.6 GeV and 38.8 GeV are displayed. We notice a 30% 
to 40% difference in normalization, which could indicate that the theoretical predictions 
calculated with fixed scales are not able to follow the energy dependence of the data. The 
use of optimized scales does not improve the agreement with data, since the predictions 
based on scale pt/3 or optimized scales give very similar results for all data with energy 
below V^ ~ 40 GeV. 

4 Higher-Energy Data and Theory 



In Fig. |]J, we compare two sets of ISR data at 63 GeV, from the R806 |2l|] and AFS [g2| 
Collaborations, with theoretical predictions. In contrast to the results for lower ener- 
gies, we observe a remarkable stability in the theoretical predictions for pt ^ 5 GeV/c 
(a;T^0.16): all our scale choices lead to similar predictions. Below px ^ 5 GeV/c, the 
predictions diverge. The choice pt/2 and the use of optimized scales give similar results, 
while the choice pt/3 leads to a cross sections which is smaller by a factor of two. Below 
pr ~ 6 GeV/c, no conclusion is possible because each data set favours a different scale 
choice. Given the statistical relevance of the R806 data, the choice pt/2 or the use of 
optimized scales appear to be most appropriate. 

The agreement of theory with both data sets is rather good for pt ^ 6 GeV/c, although 
the theoretical expectations tend to be somewhat higher than the data, a situation which 
is different from the fixed-target case. This fact may be related to the observation made 
above concerning the energy variation of the data-over-theory ratio. This can be checked, 
to an extent which is limited by the relatively large error bars, by comparing the fixed- 
scale theory to the three sets of R806 data at ^/s = 30.6, 44.8 and 62.8 GeV. We find that, 
at fixed xt, this ratio tends to decrease when the energy increases, in agreement with the 
E706 results. 



Turning to the UAl data on charged-pion production [16], similar comments as above 
can be made. The theoretical predictions at large pT values, above pT ^ 7 GeV/c, are 
relatively stable and agree well with the data, while, at lower pT values, the predictions 



start to diverge, thus bracketing the experimental points. 

5 Conclusions 

The fairest conclusion to be drawn from our studies of pion production in hadronic colli- 
sions is that the phenomenology of this process is not yet completely understood! Several 
problems can be identified. 

On the theoretical side, the main difficulty lies in the scale instability, which is signif- 
icant at low energy but disappears for y^^60 GeV, at least for large-enough pT values. 
Hopefully, this instability will be partly removed by resumming the large ln(l — z) terms 
associated with the fragmentation process, very much in the way the resummation of 
threshold factors improved the predictions of prompt-photon production at large trans- 
verse momentum [1^1^,0]. This improvement would be important in order to better 



understand the energy dependence of the cross sections, which, as we have seen, cannot 
be fully understood when using the fixed-scale approach. In the meantime, the choice of 
optimized scales, which is equivalent to setting the scales equal to pt/3 at fixed-target 
energies and to pt/2 at upper ISR energies, seems to be the most appropriate one and is 
anticipated to give results not too far from the resummed theory (c/. the case of prompt- 
photon production). 

On the phenomenological side, one needs a second generation of fragmentation func- 
tions, based on the resummed approach and also taking account of pion spectra in (se- 
lected) hadronic collisions. From the e~'~e~ data, upon which existing parametrizations are 
based, the gluon to pion fragmentation is largely undetermined, and there are very few 
experimental points to directly constrain the quark to pion fragmentation at very large z 
values. There is probably enough flexibility in the fragmentation functions to change the 
size of the predictions for hadronic collisions by 30% or so. This point is certainly worth 
a more detailed investigation. 

Comparing theory with data, the fixed-target data are found to lie systematically above 
the theory predictions, while the ISR data are compatible with or somewhat below the 
predictions and the UAl data are in perfect agreement with the predictions. This problem 
of energy dependence of the cross section at fixed xt is particularly manifest in the case 
of the E706 data. Is this an experimental problem? Is the scale dependence of the theory 
more involved than expected? Redoing the phenomenology using the resummed approach 
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will probably help to better understand this point. 

Although confused, the case of vr'' production is still much simpler than that of prompt- 
photon production, where 24 GeV and 63 GeV data are rather compatible with the 
predictions, while 30 GeV and 40 GeV data are much larger. Considering only exper- 
iments below 40 GeV and taking into account the fact that the vr*^ data taken in those 
experiments are compatible with each other, this seems to indicate that the systematic 
errors on prompt-photon production are probably underestimated. In particular, the back- 
ground subtraction necessary to obtain the prompt-photon spectrum should be carefully 
re- assessed. 

It is important for the search of the Higgs-boson decay to two photons to understand 
the production of photons and pions at LHC energies. Given the confused situation at 
lower energies, after more than twenty years of intense experimental and theoretical efforts, 
a lot of work remains to be done to achieve this goal. 
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Figure 1: Dependence on xt of experimental cross sections for inclusive prompt-photon 
production normalized to NLO predictions based on the MRST2 J^ parton densities. The 
dotted vertical lines correspond to px = 5 GeV/c for the E706 and ISR experiments. 
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Figure 2: Differential cross section Ed?a/d?p for the inclusive production of tt^ mesons 
with pT = 7 GeV and |r/| < 0.75 in the scattering of 530 GeV protons on a fixed Beryllium 
target, evaluated as a function of the initial-state factorization scale M and the fragmen- 
tation scale Mp using the MRST2 f^ parton densities. The contours of constant cross 
section are shown in the {M,Mp) plane. At each point, the renormalization scale fi is 
evaluated from Eq. ^. 
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Figure 3: Compilation of the inclusive tt^ cross sections discussed in the text as functions 
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Figure 4: Comparison of WA70 [T^ n^ data with NLO predictions for three different sets 
of parton densities and two different scale choices. The statistical and systematic errors 
are added in quadrature. 
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Figure 5: Comparison of UA6 [11] tt^ data with NLO predictions based on the MRST2 /0/ 
parton densities. Very similar results are obtained with the CTEQ4M ^ parton densities. 



16 



> 

C5 



Pi 
CO 

b 

173 



^ E706 
MRST2 




Pt [GeV] 
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Figure 9: Same as Fig. ^ but for E706 [Kj vr*^ data 
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Figure 10: Same as Fig. [^, but for scale choice Pt/3. The horizontal lines drawn at 1.2 
and 1.7 illustrate the mutual agreement of all data sets to within ±20%. Only data points 
with pt > 4.5 GeV are kept to avoid the use of too small factorization scales in the NLO 
predictions. 
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Figure 12: Comparison of R806 / pl|/ anc? ^FS' ^^/ vr*^ data with NLO predictions based on 
the MRST2 J^ parton densities. 
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Figure 13: Comparison of UAl \1(\] charged- hadron data with NLO predictions based on 
the MRST2 J^ parton densities. 
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